Pressure Effects on Magnetically-Driven Electronic Nematic States in Iron-Pnictides 
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In a magnetically driven electronic nematic state, an externally applied uniaxial strain rounds 
the nematic transition and increases the magnetic transition temperature. We study both effects 
in a simple classical model of the iron-pnictides expressed in terms of local SO(N) spins (with 
N = 3) which we solve to leading order in 1/N. The magnetic transition temperature is shown to 
increase linearly in response to an external strain while a sharp crossover, which is a remnant of the 
nematic transition, can only be identified for extremely small strain. We show that these results 
can reasonably account for recent neutron experimental data in BaFe2As2 by C. Dhital et al pQ. 



Undoped and under-doped iron-pnictides universally 
exhibit colinear antiferromagnetic ground states (C-AF) 
with ordering wavevectors (0, ±7r) or (±7r, 0) with re- 
spect to the tetragonal iron lattice. The C-AF order 
is necessarily accompanied by an orthorhombic lattice 
distortion [2j [3]. The magnetic transition temperature 
Taf an d the structural or "nematic" transition temper- 
ature TV are closely related; TV is either equal to or 
slightly greater than Taf-, TV > Taf [2 E]. The ne- 
matic distortion breaks the C4 rotational symmetry of 
the tetragonal lattice. [4] The fact that the C-AF state 
also breaks the same C4 symmetry suggests the driv- 
ing force of symmetry breaking may be the magnetism, 
itself [5 - 7J , i.e. the broken symmetry state should be 
thought of as an electronic nematic [8]. 

A number of striking experimental observations have 
been successfully interpreted in this light. Trans- 
port measurements reveal the existence of a large, 
intrinsic anisotropy in the in-plane resistivity above 
^afEJ [10]. Similar anisotropics of various physical 
quantities have been observed in scanning tunneling 
microscopy (STM) [11 j, magnetic neutron scattering [12, 
H5] , optical reflectivity measurements [14j and angle re- 
solved photoemission spectroscopy (ARPES) jT5j [16] in 
de- twinned and (even at T > TV) strained samples. 

However, the origin of the nematic state remains con- 
troversial. As these materials are all metallic, some ap- 
proaches emphasize the role of itinerant electrons, but 
the "bad metal" character of the conducting state, the 
small size of the Fermi pockets, and the relatively high 
scale of the ordering temperatures Taf and TV suggest 
that a description in terms of localized classical spins 
(or possibly orbital moments), which neglects the itiner- 
ant electrons, may be sufficient to capture the essential 
physics. (As the ordering temperatures are tuned toward 
T = 0, where quantum effects become increasingly im- 
portant and RKKY-like induced interactions become in- 
creasingly long-ranged, it certainly becomes increasingly 
problematic to ignore the effects of itineracy.) At a min- 
imum, the large size of the ordered moments (which can 
exceed l/i# at low T), and the persistently commensu- 



rate character of the ordering rules out a picture of the 
ordered state based on a weak-coupling description and 
Fermi-surface nesting. It is also open to debate the ex- 
tent to which lattice effects (e.g. electron-phonon cou- 
pling) and orbital ordering are essential drivers of the 
physics. For instance, a small but evident difference in 
occupancy of the d xz and d yz in the nematic state has 
been observed by ARPES [15] [16] in de- twinned samples 
under strain. However, it follows from symmetry that 
any correlation function which transforms like the ne- 
matic order parameter will develop a non-zero expecta- 
tion value in the nematic state, whether it is essential to 
the mechanism, or simply responding parasitically to the 
broken symmetry. 

Recent neutron scattering data from BaFe2As2 by C. 
Dhital et alpQ show that the C-AF magnetic transition 
can be affected by relatively small strain fields. In this 
paper, we adopt the most economical model which pos- 
sesses the requisite ordered phases consisting of classical, 
localized, SO(N) spins (with N = 3 corresponding to the 
physically relevant Heisenberg case) residing on the Fe 
lattice with appropriately chosen antiferromagnetic cou- 
plings. We solve this problem to leading order in 1/N in 
the large TV limit, including the effects of a small, exter- 
nally imposed uniaxial strain. We find that in response 
to a small uniform strain of magnitude Aq, the magnetic 
ordering temperature shifts [17] according to 
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where the susceptibility exponent 7 = 2 + G(l/N) , and 
X ~ (Tj^-TafY 1 as TV -> t af- So long as TV > Taf, 
the rounding of the nematic transition occurs on a scale 
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where, again for N 00, x = 1 + 0(1/ N). That we 
obtain results that satisfactorily account for the observa- 
tions of Dhital et al supports the notion that this mini- 
mal model captures much of the essential physics of mag- 
netism and nematicity in these materials. 
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Model: We start with the previously considered J\ — 
J 2 — J z — K model of magnetism in iron-pncitides: 

H = ^ ] \J\Sr,n ' Sr+6i,n ~ ^(^r,n ' Sr+5 lj n)' 2 ] 
n,f,52 n,r 

in which 5f? n is a spin S operator on the site f in plane 
n and Si and £2 are the first and second nearest neighbor 
lattice vectors in the Fe plane. J\ and J 2 are the in- 
plane nearest neighbor (NN) and next nearest neighbor 
(NNN) magnetic couplings respectively, J z is the cou- 
pling between layers along c axis and K is the NN bi- 
quadratic coupling. The C-AF groundstate arises for J 2 
sufficiently large compared to J\. (This condition re- 
duces to J 2 > Ji/2 in the limit S — >• 00.) The origin of 
K term has been discussed in [THJ [19] . 

To understand the finite temperature properties of this 
model analytically, we take the continuum limit and de- 
rive an effective classical field theory which captures the 
low energy physics of the above Hamiltonian[5 
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where 0^ is a real TV = 3 component vector field of 
unit norm [(j) n (r) ■ (j) n {r) — 1] representing the local ori- 
entation of the staggered magnetization on plane n and 
sublattice a = 1,2 as defined in ref.[5]. Here the cou- 
plings are related to those in the Hamiltonian according 
to p oc J2,t? oc Ji, p z ex J z and g oc K + AK, where 
AK ~ 0.13 Jf/SJ 2 [l + 0(1/5)] is generated by quan- 
tum fluctuations of the spin. An external strain fields 
A n (f) induces a coupling between the sublattices: 

H^H-N / d 2 r e Mf^ <f)-$? if (3) 

As a final step, we decouple the quardic term so that 
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where the Lagrangian for each plane is given by 
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where A"(r) are the Lagrange multiplierfields which en- 
force the normalization of (j) and cr n (f) is the Hubbard- 
Stratonovich fields. For this action, the nematic order is 
given as (a n ) = j^($n' > • $n^) and the magnetic order 
parameter by (<$^) = [sign((a n ))] (</>n } ). 

The layered nature of the materials is reflected in the 
fact that we will always assume p ^> p z (in units in which 
the spacing between planes is 1), and we will shortly take 
the fact that the C-AF phase is most stable for J2 ^> J\ 
to justify neglecting the effects of the gradient coupling 
between the two sublattices, we will set n = 0.[20 The 
coupling constant g determines the extent of separation 
between T/v and Taf, which empirically is small implying 
that g, too, can be considered to be small. To make this 
problem tractable, we will treat 1/N = 1/3 as a small 
parameter, i.e. we will report explicit results in the limit 
N — )> 00. Generally, N = 3 is large enough that no 
qualitative errors, and only small quantitative errors (i.e. 
in values of the critical exponents) are expected. 

For a constant strain field A(f) = A, in the TV — >> 00 
limit, the nematic order can be obtained by finding the 
saddle point of the above Lagrangian, resulting in the 
following self-consistent equations for A and a: 
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where A is momentum cutoff, and 
G- 1 = 

pk 2 + p z cos(k z ) + 2AT -(27V + 2A + 2nk x k y ) 
-(Tcr + 2A r]k x ky) pk 2 + p z cos(k z ) + 2AT 

To simplify the calculations [20 , we set 77 = 0, in which 
case the integrals can be evaluated to yield 
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where the a' = a + Aq/T, p = pA 2 /T and p 2 = p z /T. 

Solutions: In the absence of Aq, for any ^ 7^ there are 
two transition temperatures as shown in [5] . The nematic 
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transition temperature T/v is determined by the discon- 
tinuity of the function and the magnetic transition 
temperature Taf is determined by A = cr + p z . However, 
when Aq ^ 0, there is no discontinuity in the function 
because the external strain field already breaks the rota- 
tional symmetry [21]. Typical plots a and ^ as a func- 
tion of T for different values of Aq are shown in Fig(I] 
and Fig(2j When Ao is small, a crossover temperature 
Tjfr can still be identified as the inflection point temper- 
ature at which has a maximum. When Aq is large, 
^ increases smoothly as the temperature is lowered, so 
a well defined crossover temperature cannot be identi- 
fied. From the numerical data, we see that Aq ~ 10~ 3 p 
is sufficiently large to eliminate the inflection point. 

FIG. 1. Plot of the nematic order parameter (a) vs. 
temperature for different values of external field Ao — 
(red), 0.001 (green), 0.01 (blue). The parameters chosen are 



pA = 1, g = 0.3 and p z 



0.01. 
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FIG. 2. Plot of the derivative of the nematic order parameter 



vs. temperature for different external fields Ao(10 
parameters used are p z — 0.01, g — 0.3, p — A = 1. 

do- 
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Taking A = <j a +p z /T, we can derive the magnetic tran- 



sition Taf, which is determined by the following equa- 
tions, 
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p z , the magnetic transition temperature T\ F is approxi- 
mately 
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where T^ F 



For small external strain field, 
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ture, AT af =T A f(Ao) 
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Since oaf ~ <?, the coefficient in the right side of the 
above equation goes as ~ 1/g for small g. The above 
results can be extended even in the limit of g — >> 0. It 
is easy to show that in this limit, T/v — Taf ~ # 1//2 - 
Plugging this into EqjHJ we obtain the expression below 
Eq-B 

Eq. [I] can be further checked in the case of g = 0. For 
g = 0, the change of magnetic transition temperature is 
given by 
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consistent with the expectations of scaling theory, Eq. [T] 
Comparison with experiment: Fig. 3 shows the com- 
parison of our theoretical results with experimental ob- 
servations of Dhital et al for the magnetic transition and 
nematic crossover as a function of uniaxial strain. The 
parameters used to generate the theoretical results, rep- 
resented by the lines in the figure, are presented in the fig- 
ure caption. The upper curve, representing the nematic 
crossover, is solid where there is a well-defined inflection 
point associated with the crossover, and a dashed line 
where the crossover has become so smooth that there 
is no local maximum in the temperature derivative of 
the nematic order parameter. The lower line indicates 
the theoretical magnetic ordering temperature. Fig. 4 
shows the variation of the nematic transition with Ao 
from which Fig. 3 was obtained. The close correspon- 
dence between the theoretical and experimental curves 
supports the conjecture that the starting model captures 
the essential physics, although the comparison involves 
too many empirically determined parameters to make 
this conclusion inescapable. In order to match the ex- 
perimental transition temperature, the value of p used in 
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FIG. 3. Comparison of the external field, Ao, variation 
of Magnetic Ordering Temperature and Nematic transition 
temperature with experiment. The parameters chosen are p z — 
0.06p,# = 0.3p,p~71if. 
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FIG. 4. Variation of the Nematic order jump with external 
field Ao. The parameters chosen are the same as Fig. 3 with 
N=3 and A in units of 10~ 3 . 




our theoretical calculation ( see the caption of Fig. 3) 
is smaller than the value measured in neutron scattering 
experiments [12 . This is expected since the transition 
temperature is always overestimated when N is taken to 
be 3 in the large N expression. 

In summary, we have shown that the minimal model 
with short ranged magnetic exchange couplings can sat- 
isfactorily account for the change of both structural and 
AF transition temperatures under the uniaxial stain mea- 
sured by Dhital et al in BaFe2As2. The experimental 
results strongly support the notion of the magnetically- 
driven nematicity in iron-pnictides. 
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